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I.  PROBLEM  STUDIED  AND  RESULTS  OBTAINED  UNDER  THE  ARMY  RESEARCH  GRANT 


The  research  supported  by  the  Army  was  to  study  some  aspects  of 
three-dimensional  waves  in  inelastic  media.  Inelastic  media  include  elastic- 
plastic  materials  and  viscoelastic  solids.  The  following  is  a  summary  of 
the  results  obtained. 

Basic  information  needed  in  analyzing  wave  propagation  in  three- 
dimensional  media  are  the  eigenvectors  (or  the  polarization  vectors)  of  the 
acoustic  tensor  and  the  transport  equations  which  determine  the  growth  or 
decay  of  a  discontinuity.  These  were  obtained  for  a  fairly  general 
anisotropic  elastic-plastic  media  in  which  the  eigenvectors  in  the  plastic 
state  were  expressed  explicitly  in  terms  of  the  elastic  acoustic  tensor  [1]. 
This  expression  makes  it  easier  to  calculate  numerically  the  eigenvectors 
in  the  plastic  state.  The  discontinuity  in  acceleration  aoross  the  singular 
surface  is  proportional  to  this  eigenvector  regardless  of  whether  the 
singular  surface  is  an  elastic-elastic,  plastic-plastic  or  elastic-plastic 
interface.  We  also  obtained  the  ray  velocity  and  the  transport  equations 
for  the  growth  or  decay  of  the  discontinuity  along  the  ray.  The  possibility 
of  shock  waves  in  a  plastic  region  is  explored.  Finally,  the  general 
results  are  specialized  for  isotropic  elastic-plastic  solids. 

With  the  general  results  obtained  in  [1],  we  then  study  basic  problems 
needed  for  three-dimensional  waves.  The  first  is  the  wave  speeds.  Although 
the  plastic  wave  speed  in  one-dimensional  media  is  well-known,  the  plastic 
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the  direction  of  propagation,  even  when  the  material  is  isotropically  work¬ 
hardening.  However,  regardless  of  the  direction  of  propagation,  it  was 
found  that  the  difference  between  c^,  the  fast  plastic  wave  speed,  and 
Cj,  the  fastest  elastic  wave  speed,  is  no  more  than  3S%  of  c^  for  the 
von  Mises'  materials  and  29%  for  the  Tresca's  materials  [2,3].  For  materials 
which  are  not  ideally  plastic  and  for  loading  of  a  magnitude  such  that  the 
tangent  modulus  at  the  stress  state  is  not  nearly  zero,  the  difference  is 
smaller.  The  difference  will  be  even  smaller  if  the  direction  of  propagation 
is  not  in  the  direction  of  the  minimum  wave  speed.  In  fact,  one  can  always 
choose  a  family  of  directions  such  that  c^  and  c^  are  identical.  These 
results  apply  also  to  ideally  plastic  solids. 

The  significance  of  this  finding  is  that  when  we  analyze  wave  propaga¬ 
tion  in  a  thin  rod  as  a  three-dimensional  circular  cylinder,  the  difference 
between  c^  and  the  statistical  average  of  c^  for  all  directions  will 
probably  be  within  10%  of  c^.  Since  the  precursor  wave  front  propagates 
at  the  fastest  wave  speed,  this  implies  that  the  plastic  wave  front  will 
propagate  at  a  speed  within  10%  of  the  elastic  wave  front.  Thus  the 
phenomenon  observed  in  the  experiments  that  the  plastic  wave  front  propagated 
at  the  elastic  wave  speed  might  not  be  entirely  due  to  the  rate  effects  as 
widely  believed.  The  approximation  of  three-dimensional  circular  cylinder 
by  a  one-dimensional  rod  may  also  be  responsible,  if  only  partially,  for  the 
phenomenon. 

The  second  basic  problem  is  the  reflection  and  transmission  of  a  wave 
front  at  an  interface  between  two  different  materials.  A  special  case  of 
this  problem  is  the  reflection  of  a  wave  front  from  a  boundary  which  may  be 
fixed  or  stress  free  [4].  Near  the  point  of  reflection,  the  geometry  of  the 
boundary  and  the  incident  w?  e  front  is  not  important  as  far  as  the 


acceleration  waves  are  concerned.  Therefore,  we  analyzed  the  reflection  of 


acceleration  waves  from  the  surface  of  a  semi-infinite  elastic-plastic  medium 


Depending  on  the  nature  of  the  incident  waves  and  the  incident  angle,  there 
may  be  one  to  four  reflected  waves  generated.  The  analyses  can  predict  the 


number  and  the  speeds  of  waves  reflected  and  also  identify  whether  the 
regions  behind  the  reflected  waves  are  either  elastic  or  plastic.  This 


information  serves  useful  guides  for  solving  numerically  three-dimensional 


wave  propagation  problems 


Another  basic  problem  for  wave  propagation  in  elastic-plastic  media 


is  the  question  of  discrepancies  between  the  rate- independent  theory  and 


the  experimental  observations.  The  rate-effect  has  been  considered  as  the 


main  cause  of  the  discrepancies  but  there  are  also  other  factors,  such  as 


the  three-dimensional  effect.  To  study  this  effect,  we  consider  a  three 


dimensional  numerical  solution  to  wave  propagation  in  a  semi-infinite 


circular  cylindrical  rod  of  elastic-plastic  rate- independent  material  due 


to  a  uniformly  distributed  axial  load  applied  at  the  end  of  the  rod  [5] 


The  experiments  showed  that  the  wave  front  propagated  at  the  speed  of 


elastic  bar  velocity  c  while  the  rate- independent  one-dimensional  theory 


predicts  that  the  disturbance  propagates  at  the  plastic  wave  speed  c 


The  three-dimensional  numerical  solution  reveals  that  while  the  stress 


history  at  a  station  10  radii  from  the  impact  end  agrees  qualitatively 
with  the  one-dimensional  theory,  the  wave  front  travels  at  a  speed  some 


where  between  c  and  c  .  In  view  of  the  oscillatory  nature  of  the 


numerical  solution,  it  is  difficult  to  pinpoint  the  exact  arrival  time  of 
the  wave  front.  However,  the  wave  front  definitely  travels  faster  than  c 
In  [1],  we  discussed  the  possibility  of  generating  a  shock  wave  in 
three-dimensional  elastic-plastic  media.  Before  we  proceeded  to  study  in 
more  detail  the  nature  of  three-dimensional  shock  waves,  we  looked  at  the 
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existing  publications  on  one-dimensional  shock  waves  and  discovered  some 
interesting  results.  Firstly,  while  the  infinitesimal  deformation  theories 
are  adequate  if  the  strains  are  small,  they  lead  to  physically  unacceptable 
phenomena  if  the  impact  velocity  is  transonic  or  supersonic  [6] .  Secondly, 
one  may  define  the  growth  of  a  shock  wave  by  the  growth  in  the  discontinuity 
in  the  velocity  across  the  shock  wave  as  the  shock  wave  propagates.  One 
may  also  define  the  growth  of  a  shock  wave  by  the  growth  in  the  discontinu¬ 
ity  in  the  stress,  strain  or  the  entropy.  It  turns  out  that  while  one 
definition  predicts  the  growth  of  the  shock  wave,  the  other  definition  may 
predict  its  decay  [7] . 

The  problem  of  three-dimensional  shock  waves  in  solids  appears  to  be 
much  more  complicated  than  we  had  anticipated.  We  were  not  able  to 
accomplish  the  study  during  the  period  of  the  grant.  However,  we  were  able 
to  re-examinq  the  three-dimensional  shock  waves  in  elastic  fluids  [8] . 

Knowing  that  we  could  not  complete  the  study  on  three-dimensional 
shock  waves  in  solids,  we  studied  the  problem  of  wave  propagation  normal 
to  the  layering  of  a  periodically  layered  elastic  or  viscoelastic  medium. 
Even  though  many  approximate  theories  have  been  proposed  for  wave  propaga¬ 
tion  in  a  composite,  there  appears  to  be  no  reliable  way  to  predict 
satisfactorily  the  transient  response  in  the  region  which  is  neither  far 
away  from  the  impact  end  nor  near  the  wave  front.  The  same  is  true  for 
transient  waves  in  a  finite  layered  medium.  We  proposed  a  theory  which, 
from  a  practical  point  of  view,  can  determine  satisfactorily  the  transient 
solution  in  a  region  which  can  be  near  or  far  away  from  the  impact  end  [9] . 
The  fundamental  idea  of  the  theory  is  to  replace  the  layered  medium  by  an 
"equivalent"  homogeneous,  linear  viscoelastic  medium  whose  dynamic 
response  is  identical  to  that  of  the  layered  medium  at  points  which  are  the 


centers  of  the  odd  layers.  Since  wave  propagation  in  a  homogeneous,  linear 
viscoelastic  medium  can  be  solved  easily  by  many  known  numerical  schemes, 
one  can  solve  any  transient  wave  propagation  problem  in  semi-infinite  or 
finite  layered  medium  [10]  by  the  present  theory.  Several  numerical 
examples  are  obtained  and  compared  with  the  exact  solutions  by  the  ray 
theoy.  The  agreement  between  the  two  solutions  is  extremely  satisfactory. 
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